The dependence of the dry friction threshold on rupture dynamics 
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The static friction coefficient between two materials is considered to be a material constant. 
We present experiments demonstrating that the ratio of shear to normal force needed to move 
contacting blocks can, instead, vary systematically with controllable changes in the external loading 
configuration. Large variations in both the friction coefficient and consequent stress drop are tightly 
linked to changes in the rupture dynamics of the rough interface separating the two blocks. 

PACS numbers: 46.55.+d, 46.50.+a, 62.20.Qp, 81.40.Pq 



For centuries the onset of frictional motion has been 
described by the concept of a friction coefficient, /i, that 
reflects proportionality between applied shear and nor- 
mal forces when frictional motion initiates. This concept, 
which dates back to DaVinci, Amontons, and Coulomb, 
was first explained by Bowden and Tabor [l| as resulting 
from the proportionality of the real area of contact, A 
with the normal force, Fn, pressing two bodies together. 
Phenomenological corrections describing the dependence 
of the friction coefficient on time and sliding velocity were 
later introduced [2|-l5|. Although significant, these correc- 
tions generally predict variations of /i of a few percent. 

The onset of frictional slip is, however, caused by the 
rapid rupture of the discrete ensemble of contacts that 
forms a rough frictional interface The nature of 

these rupture modes is coupled to the local profiles of 
the shear, r(x), and normal, cr(x), stresses along the in- 
terface with local variations in r(x)/a(x) significantly in- 
fluencing the rupture mode selection [8[. Moreover, for 
given external loading conditions r(x)/a(x) can not only 
vary considerably but can, locally, surpass macroscopic 
friction coefficient, /i, by hundreds of percent [8[, demon- 
strating that ji is not the local criterion for slip onset. 

We will demonstrate that the static friction coefficient, 
/i5, is not a constant material property. It, in fact, varies 
systematically by nearly a factor of two with the exter- 
nal loading configuration. Furthermore, the threshold 
for frictional sliding is intimately linked to the rupture 
dynamics of the interfacial contacts. This link provides 
a mechanical explanation for significant variations of fis 
found even within a single experiment with a set external 
loading configuration. 

The experiments described here were all in the stick- 
slip regime, where each experiment consisted of a se- 
ries of stick-slip events. The experimental system, de- 
picted schematically in Fig. [TJ is similar to that de- 
scribed in HQ. Two PMMA blocks (shear wave speed, 
Cs = 1370m/ s) were pressed together with a uniform 
normal force Fn. The top (bottom) block dimensions 
were 200 x 6 x 100mm (300 x 60 x 30mm), where x, y and 
z are, respectively, the propagation, thickness and nor- 
mal loading directions. The bottom block was mounted 
on a low-friction translational stage constrained to move 
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FIG. 1. (a) Schematic view of the experimental setup. Nor- 
mal force Fn is applied uniformly to the top block. Shear 
force Fs is applied to the translational stage parallel to the x 
axis. Stresses were measured along the frictional interface by 
miniature strain-gages (green squares, not to scale). A small 
tilt, a, was introduced to controllably vary the stress profiles 
(call-out), (b-left) The initial cr(x) profiles for a = 0.01° and 
Fn = 2000 — 7000iV (see legend), showing how stress profiles 
differ with Fn. (right) a(x) profiles immediately prior to slip 
initiation for events from experiments with these values of Fn 
(color coding as in (b)). For each Fn, the initial a(x) profiles 
are changed to a different degree by the Fs induced torque. 



along the x axis. Relative changes in the real contact area 
along the rough interface, A(x,t), during each stick-slip 
event were measured by a method based on total inter- 
nal reflection 0, E| Q , where t denotes the time following 
the rapid rupture initiation for each event. A(x,t) are 
normalized by values lms prior to rupture initiation at 
t = 0. 

Shear force, Fs, was applied to the bottom block 
in the negative x direction via a load-cell of stiffness 
10 7 A7m at a loading rate of 10 - 20N/s. F N and F s 
were continuously measured at 100Hz. In addition, Fs 
was measured at 250KHz during stick-slip events to ac- 
curately determine its peak value for each event (each 
event corresponded to a sharp drop in Fs). Force ra- 
tios fi(t) = Fs (t) I Fn (t) , were measured to better than 
1% accuracy. For each event, we define [is = fi(t = 0). 
Rapid acquisition of A(x, t) and Fs was triggered by an 
acoustic sensor mounted to the x = face of the top 
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block. Local strain values were measured ~ 2mm from 
the interface [8[, yielding profiles of a(x) and r(x) prior to 
each slip event, at 56 and 8 locations along the interface, 
respectively. 

To controllably vary both the stress profiles along the 
interface and the resulting rupture dynamics, we varied 
the loading conditions by slightly tilting the top block 
at angles a = 0.01° and a = 0.02° (see Fig. ffl[a)). 
Note that while the absolute value of a is accurate to 
within 10%, its value was reproducible to within 1% in 
sequential experiments. For a full discussion of loading 
see jsj. For a given a, variation of Fn over the range 
2000 — 70007V controllably changes the stress profiles. 
This is illustrated in Fig. [ljb-left) which shows how <j(x) 
changes as Fn is increased for a = 0.01°. With a in this 
orientation, these changes in <j{x) are partially offset by 
the torque applied by Fs- Fig. QJb-right) shows how Fs 
influences <j(x) profiles at t — 0. a(x) profiles vary signif- 
icantly for different Fn at t = due to different degrees 
of compensation by the F^-induced torque. Increasing 
a changes both the initial profiles and resulting stress 
profiles at slip initiation. 
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FIG. 2. (a-c) Measurements at 20/is intervals of the changes 
in the real contact-area A(x,t) for representative events 
with a = 0.01° and (top-bottom) F N = 2000,4000 and 
6000iV. Hotter (colder) colors indicate increased (decreased) 
A. A(x,t) was normalized by its value 1ms before rupture 
initiation. Each rupture consists of a slow front initiating 
near x = which sharply transitions, at x = It to rapid 
propagation. Transitions to front velocities V = 2000, 1800 
and 1500m/s occur at It = 27,38 and 52mm for (a-c). 
/is decreases with increasing It. (d) /is vs. It normalized 
by the interface length L — 200mm, for 122 events with 
a — 0.01° (open symbols) and a — 0.02° (filled symbols) 
for Fn = 2000 — 7000 N. For each a, all points collapse to 
distinct well-defined curves, each offset in /is- Whereas data 
for a = 0.01°, values of /is roughly cluster for each Fn, data 
for both F N = 7000iV and a = 0.02° span most of the col- 
lapse curves, indicating that /is depends on rupture dynamics 
rather than directly on Fn- 

In accepted views of static friction, [is is a material 



property which should not depend on Fn- Fig. [2] demon- 
strates that /is, in fact, changes significantly with the 
loading details. Furthermore, the variation of /is is sys- 
tematically correlated with the dynamics of the interface 
rupture. Examples of how rupture dynamics are affected 
by changing Fn for initial values of a, are presented in 
Fig. [2j The left panel presents A(x,t) measurements 
for three representative slip events. The front dynam- 
ics for all events are qualitatively similar; a slow front 
(30 < V < 300m/s) initiates near x = and undergoes a 
sharp transition to rapid rupture (1300 < V < 2000ra/,s) 
which traverses the rest of the interface. These dynamics 
are consistent with the "nucleation phase" of [7] . Quanti- 
tatively, however, the length of each slow front at transi- 
tion, denoted by It, varies with the changing stress pro- 
files for different Fn- The resultant {is systematically 
decrease as It grows (Fig. [2fa-c)). In Fig. [2^d) [is 
versus It (normalized by the interface length, L) is pre- 
sented, for both a = 0.01° and a = 0.02°, over a range 
of Fn- Included are 122 events from 11 experiments. 

The data for a = 0.01° roughly cluster around a dif- 
ferent value of [is for each value of Fn- Us, however, 
does not explicitly depend on Fn- A close look at [is for 
both Fn = 70007V and a = 0.02°, reveals that these data 
do not cluster for specific values of Fn, but instead span 
nearly the entire range of [is- When plotted as a function 
of It, the "scattered" data from all experiments collapse 
onto well-defined curves with very little scatter for each 
a. While the two data sets are offset in [is, events with 
the same a have similar slopes. Thus, the systematic 
variation of [is is not an explicit dependence on Fn, but 
is, instead, closely tied to the rupture dynamics. 

For these loading conditions, the rupture dynamics of 
the first event in each experiment generally differ from 
those of subsequent events. Stresses along the interface 
are induced by Poisson expansion that is frustrated due 
to frictional pinning 0, E3]. These stresses often lead, 
in this configuration, to rupture nucleation of the first 
event away from the sample edges (c.f. Fig. HJ. These 
data were, therefore, not included in Fig. [2]since, in these 
cases, It is not defined. 

Large variations in both the peak values of [± ([is) and 
the consequent drops, A/i, in each stick-slip event are ap- 
parent in the representative time series of successive slip 
events presented in Figl^a). [is variations are notably 
larger (< 8%) for a = 0.02° than for a = 0.01° (< 3%) 
for the same Fn. The size of these variations are not 
explained, for example, by logarithmic aging in PMMA, 
f^s(t) = Ms'(O) + 0.01log(t) [3|, |4(, which accounts for only 
a ~ 0.4% variation in /is over the temporal-interval vari- 
ations between slip events. In general, distributions of 
both [is and A/if/is are quite broad, as demonstrated in 
FigEl This apparent random aspect of frictional stabil- 
ity disappears when both [is and Afi/fis are plotted as 
a function of It- In particular, Fig. [3fc) demonstrates 
that the broad distribution of Afi/fis (Fig. [3fc-inset)), 
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obtained for all experiments with both a values, collapses 
onto a relatively well-defined curve, when plotted as a 
function of It- Thus, the aperiodic nature of stick-slip 
friction is also closely related to rupture dynamics. 
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FIG. 3. (a) The loading curve a as a function of time for 
F N = 300(W and a = 0.02°. /i(t) drops by A/i after each 
event. Values of /is (triangles) vary by up to 8% in subsequent 
events, /is and the relative stress drops, A/i//is, collapse to 
well-defined curves when plotted vs. It/L. (b) /is vs. It/L 
(inset) A/i/fis vs. It/L for the data in (a), (c) A/i//is vs. 
It/L for the broadly distributed Aji/jis (inset). Plotted are 
all data from a = 0.01° - 0.02° presented in Fig. [2]i. 

To this point we have seen how controlled changes in 
the loading configuration (either different Fn and/or a) 
can, by varying the rupture dynamics, induce system- 
atic changes in /is of order 10%. To demonstrate the 
dependence of /is on rupture dynamics, we used partic- 
ular examples of loading configurations in which rupture 
dynamics could be roughly characterized by a single pa- 
rameter (e.g. It). 

Even in these simple cases, It is not uniquely deter- 
mined by details of the loading (c.f. Fig. 2). In the 
general case, moreover, the stress profiles and ensuing 
rupture dynamics are strongly dependent on the loading 
configuration [8] and are not, necessarily, characterized 
by a single parameter. In Fig. [4j we show that the spread 
in /is can, in fact, be much larger when a larger variation 
of stress profiles along the interface is imposed. Fig. 2] 
presents /is data for 458 slip events that resulted from a 
variety of different loading configurations. These range 
from pure edge loading (where Fs is applied directly to 
a sample edge) to the loading configurations described in 
Figs. HEl Further details of the various configurations 
can be found in [4|, [8j. The figure shows that, simply by 
changing the loading configuration, /is can vary by a fac- 
tor of almost two for the same two blocks. Note that Fig. 
Ufa) is not a histogram, but simply a record of the num- 
ber of events in which a given /is was measured. Sub- 
populations of events from single experiments (colored 
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FIG. 4. (a) Count of events with varying /is, from exper- 
iments conducted under widely varying loading conditions. 
These range from pure edge loading to near uniform appli- 
cation of Fs. Note that this plot is not a histogram and 
solely reflects the number of experiments performed whose 
loading conditions yielded /is for each range. The color bars 
(bar width = two standard deviations) in the plot depict sub- 
populations of events from single experiments, each with dif- 
ferent loading conditions. The bar at /is = 0.62 corresponds 
to Fig. 3(a-b). (b-c) A(x,t) (top) and stress ratio profile im- 
mediately prior to the event (bottom) for two experiments at 
extreme values of /is- (b) Applying Fs to the edge of the 
sample resulted in strong localization of r near the edge, and 
initiation of a slow (~ 90m /s on average) rupture which tra- 
versed most of the sample, resulting in /is = 0.4. (c) The 
relatively uniform stress-ratio profile with lower stress-ratio 
near the edges inhibited rupture initiation, resulting in a su- 
pershear rupture initiating at x ~ 158mm with /is = 0.73. 



bars in Fig. 11(a)) demonstrate that for given loading 
configurations, /is is relatively well-defined. For given 
loading conditions, the maximal spread of /is is close to 
the ~ 8% variation presented in Fig. [3j 

In Fig. |4^b),(c) we present two sample slip events 
which correspond to the low and high extremes of /is 
within Fig. @|a). The /is = 0.4 event (Fig. 0(b)) was 
obtained by application of Fs solely to the edge of the 
top block. This lead to a r/a profile whose maximal 
value was highly localized near x = 0, and much lower 
elsewhere. In [8] it was shown that once a rupture front 
nucleates, ruptures will propagate in any of three differ- 
ent modes, as long as r(x)/a(x) locally exceeds ~ 0.25. 
Here, a rapid front (V ~ OACs) nucleated at the high 
r/a ratio at x = and then continued at x « 60mm as 
a slow rupture front (V ~ 0.07Cs) across most of the in- 
terface, before transitioning back to rapid propagation as 
the rupture approached the x = 200mm edge. As most 
of the interface ruptured near the minimal values of r/a 
(~ 0.25) that support rupture propagation, the overall 
value of /is, which is the ratio of the spatial integrals of 
r and a, was extremely low. 
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In the us = 0.73 event (Fig. 13(c) ) the values of 
r(x)/a(x) across nearly the entire interface were both 
significantly higher (~ 1) than in the jus = 0.4 event 
and uniform, except for the relatively low values near 
the edges. In this case, the rupture nucleated within 
the interior of the interface, at a point (x = 158mm) 
where r/a was near maximal. Once nucleated, the result- 
ing counter-propagating ruptures sharply transitioned to 
super-shear propagation (V ~ 1.85Cs), as predicted in 
[8[. Here, the high overall value of /is reflects the large 
mean values of the r/cr spatial profile. A cardinal differ- 
ence between the two events in Figs. 13(b), (c) is in their 
rupture nucleation threshold. 

We have demonstrated that global quantities (e.g. fig 
and stress drops) are both strongly dependent on the im- 
posed loading configurations and tightly linked to rup- 
ture dynamics. We suggest that the nucleation thresh- 
old for a given loading configuration may be analogous 
to the onset of fracture. Initiation of rapid fracture (the 
Griffith criterion) will occur when the amount of bulk po- 
tential energy released by extension of an existing crack 
surpasses the fracture energy, the energy required to cre- 



ate the new surfaces needed to extend the crack [12 1. 
Fracture initiation, in any given loading configuration, is 
solely determined by this energy balance and not by the 
force balance inherent in the concept of a "friction co- 
efficient". In many loading configurations, the fracture 
threshold strongly decreases with the increased length of 
an initially imposed "seed" crack. The systematic de- 
crease in jiis with It in Figs. EE] could suggest that the 
length of the slow front phase (It) plays an analogous 
role to this initial seed crack in determining the stabil- 
ity of a frictional interface to rapid rupture. The role 
of r/a in this process is also understood in this analogy. 
a(x) oc A(x) is related to the local frictional strength, 
analogous to the fracture energy. In contrast to fracture, 
the frictional strength is not a material constant; it varies 
spatially and temporally with A(x,t). r(x) reflects the 
effective amount of stored potential energy surrounding 
a rupture tip. 

The above analogy provides a conceptual framework 
for understanding how and why the onset of frictional 
slip can vary. It, furthermore, relates {is and rupture 
dynamics. It is, however, far from the entire story. Un- 
like fracture, where dissipative processes only occur at a 
crack's tip and the crack faces are stress-free, in frictional 
slip the "faces" behind the rupture tip are far from free. 
Although the frictional strength in the wake of the rup- 
ture tip is reduced, nearly the same level (80%) of contact 
area persists behind the rupture tip 0, El [l3| during the 
rupture process. A theory that fully accounts for these 
important effects does not yet exist. 

This intuitive framework also provides a "mechanical" 
explanation for widely different measured values of /is 
obtained for the same materials [H, 0, S wnen dif- 



ferent loading conditions are applied. This additionally 
explains differences in the local and global friction co- 
efficients observed in simplified models [Hj. Dynamical 
changes to the stress profiles along an interface in the 
course of a single experiment also may explain the origin 
of aperiodic stick-slip behavior fllj . 

We are not yet able to predict the onset of slip for a 
given configuration. Our results, nonetheless, indicate 
that the answer may well lie in obtaining a better under- 
standing of the evolution of the stress profiles driving the 
frictional rupture. This is a fundamental open question; 
of paramount importance to both mechanical stability 
and earthquake prediction. 
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